Topological spin liquid on the hyper-kagome lattice of Na^^Os 
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Recent experiments on the "hyper-kagome" lattice system Na^rsOs have demonstrated that it is 
a rare example of a three dimensional spin- 1/2 frustrated antiferromagnet. We investigate the role 
of quantum fluctuations as the primary mechanism lifting the macroscopic degeneracy inherited by 
classical spins on this lattice. In the semi-classical limit we predict, based on large- N calculations, 
that an unusual q = coplaner magnetically ordered ground state is stabilized with no local "weather 
vane" modes. This phase melts in the quantum limit and a gapped topological Z2 spin liquid phase 
emerges. In the vicinity of this quantum phase transition, we study the dynamic spin structure 
factor and comment on the relevance of our results for future neutron scattering experiments. 



Introduction. Geometrically frustrated spin systems 
are a promising place to search for new "exotic" phases 
of matter. Classical spin models on such lattices typi- 
cally exhibit massive degeneracy of their zero tempera- 
ture ground state ensembles so that phases of matter in 
these systems are selected exclusively by thermal and/or 
quantum fluctuations. While magnetically ordered states 
may still arise at very low temperatures in such systems, 
emergent phases with, for example, fractionalized excita- 
tions and topologically protected ground states, are also 
possible, especially in spin ^ systems p], Q. 

On the other hand, spin 1/2 frustrated magnets are 
extremely rare, with only a few known examples such 
as Volborthite[3], Herbertsmithite[4] and certain organic 
salts [El, In light of this, the recent experiment Q 
demonstrating that Na^raOg (NIO) is a new spin 1/2 ge- 
ometrically frustrated magnet and the first such material 
to contain a genuinely three-dimensional lattice is espe- 
cially important. In particular, these experiments see no 
sign of magnetic or orbital ordering down to a few Kelvin 
while bulk susceptibility measurements show that the 
spin 1/2 iridium atoms interact anti-ferromagnetically 
with a Curie- Weiss constant of 650 K. Furthermore a 
nearly constant susceptibility and a C/T that is mag- 
netic field independent with a broad anomalous peak at 
about 25K is observed. Certainly, the remarkable hyper- 
kagome lattice (see Fig. [1]) that spin 1/2 iridium atoms 
occupy is responsible for these striking observations. 

In order to investigate whether NIO's magnetic prop- 
erties are truly quantum mechanical in origin, it is nec- 
essary to understand how the classical ground state de- 
generacy is lifted by both thermal fluctuations (entropic 
selection) and quantum fluctuations (energetic selection). 
Fortunately, much is now known about the entropic se- 
lection of these classical ground states due to a recent 
study of the classical Heisenberg model on NIO's hyper- 
kagome lattice [8]. In particular, the classical cooperative 
paramagnetic (CCP) phase exists over a wide temper- 
ature range with a spin nematic phase setting in be- 
low J/1000 w 300 mK, where J is the Heisenberg ex- 
change coupling. However, the behavior of the classical 
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FIG. 1: (a) The unit cell of the iridium atoms in NIO. All sites 
and bonds are chemically equivalent due to the P4i32 space 
group symmetry of the lattice. Letters indicate the 120° co- 
planer magnetic ordering expected in the semi-classical limit, 
(b) Positions of magnetic Bragg peaks in the ZiO^-plane of the 
ordered phase shown in (a). Notice magnetic Bragg peaks for 
this phase coincide with some lattice Bragg peaks. 



model deviates from measurements of NIO's thermody- 
namic quantities at and below the anomalous peak in 
C/T at 25 K suggesting that quantum fluctuations may 
dominate even at relatively high temperatures. 

In this letter, we present litmus tests for the ex- 
istence of a quantum spin liquid phase in NIO. We 
do so through an investigation of the ground states of 
the Heisenberg quantum antiferromagnetic spin model 
on its frustrated hyper-kagome lattice with the large-TV 
Sp(7V) method[l, QjJ, [Tl|. This allows us to study the 
semi-classical (large "spin" ) and quantum (small "spin" ) 
regimes on equal footing. In particular, we show how a 
quantum spin liquid phase in the quantum regime can be 
distinguished from both a CCP phase at finite tempera- 
tures and the most likely magnetically ordered phase. 

In the semi-classical regime, we predict a ground state 
with the 120° co-planer magnetic ordering shown in Fig. 
QJi. Defining ABC as three unit spin vectors whose 
sum vanishes, this ordering is best characterized by a 
BCBC pattern along the "threads" , paths shown by the 



darker bonds. This is in contrast to the x y/3 state 
on the two-dimensional kagome lattice which is also co- 
planer but with a BCBC pattern around hexagons. Since 
threads are as long as the length of the syst em, this state 
is not expected to be entropically selected[12], a result 
consistent with the thermally stabilized spin nematic or- 
dering and an absence of magnetic ordering found in the 
Monte Carlo study of Ref. Fig. [TJd shows the position 
of the corresponding magnetic Bragg peaks. 

On the other hand, in the quantum regime the mag- 
netic ordering melts and a Z2 quantum spin liquid phase 
arises. Since the spin liquid is a paramagnetic phase, it is 
important to distinguish it from the CCP phase. To this 
end, consider the dynamic spin structure factors S(k,uu) 
depicted in Fig. O The energy integrated spin structure 
factor of the Z2 spin liquid ground state in Fig. [2^ can 
be compared directly with the static spin structure factor 
of the CCP phase in Fig. [2]o[8]. In the CCP phase, long 
range dipolar correlations lead to the vanishing of spin 
correlations (nodes) along the [hhh] and symmetry re- 
lated directions. In contrast, dipolar correlations in the 
Z2 spin liquid phase retain a finite correlation length, 
leading to the absence of this nodal structure. 

Another important signature of the quantum spin liq- 
uid phase is the presence of spin- 1/2 neutral excitations, 
dubbed spinons. This leads to a spinon-antispinon con- 
tinuum in the dynamic spin structure factor and the bot- 
tom (or threshold in energy for a given momentum) of 
this continuum should have a well defined dispersion re- 
lation. Conversely, in a CCP phase no such sharp thresh- 
old is possible making its observation away from Bragg 
peaks a definitive proof for the existence of a quantum 
spin liquid phase. This threshold for the Z2 quantum 
spin liquid phase proposed here is plotted in Fig. [2t (the 
shaded region simply represents where S(k,cu) 7^ 0). 

The model. Consider now the Heisenberg model on 
NIO's hyper-kagome lattice. This model in the Schwinger 
boson representation is as follows [13j: 
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where the spin operators are Si = ^bi a T a ^>bi a > , the sin- 
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glet creation operator is A\j — 

constraint bl a b ia = rii = 2S is imposed on each 
site, where S is the size of the spin. It is convenient 
to generalize this SU(2) model by introducing N flavors 
of Schwinger bosons, b^ : m G {!,..., TV} and letting 
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create a singlet for each flavor. Then 
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FIG. 2: Distinguishing the Z2 quantum spin liquid ground 
state from the classical cooperative paramagnet (CCP) via 
the dynamic spin structure factor S(k,u). (a) and (b): the 
absence of the dipolar correlations in (b) lead to finite spin 
correlations along the diagonal [hhh] direction in (a) (axes run 
from h = —2 ... 2, / = —2 ... 2). (c) the spinon-antispinon 
continuum in S(k,cu): shaded (blue) region represents where 
the structure factor is finite. Notice that the threshold of the 
spinon-antispinon continuum has a well defined dispersion. 



and is Sp(7V) symmetric. In the large- TV limit Qij = 

Ylm=i(^ij) ^ s a mean fi e ld and the chemical potential 
Xi = A enforces the constraint on average where n = 
2S/N is held fixed as TV -> 00. 

This large- N limit is well suited to our purposes since 
as a function of k we can access both a quantum regime 
at small k and a semi-classical regime at large k. 

Quantum regime. To study this regime, it is useful 
to start from an extreme quantum limit and expand the 
ground state energy in powers of k. Following Ref. [3], 
before expanding in powers of ft, consider first reformu- 
lating the ground state energy in matrix form 
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where (Q)ij = Qij and then re-scaling Qij — > aQij and 
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Esp(N) scale invariant. Solving the constraint equation, 
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<9£?s p (jv)/<9A = 0, order by order in k then leads to 
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where P n = Tr[Q Q^] n /iV s ^ es is a measure of the total 
field strength, or flux (defined below), through all loops of 
length 2n starting and ending at the same site (including 
smaller loops whose total length is 2n). Through any 
individual loop C this (normalized) flux is defined by 

2 cos ($ c ) = Qij (-Q* k ) . . . Q mn (-Q* ni ) + c.c. (5) 

Thus, P n is the total flux through all loops of length 2n. 

It follows from Eq. (jlj) that the ground state must 
satisfy a flux expulsion principle: in this state the flux 
&c through all loops must vanish, if possible. To see 
that this is the case, consider using Eq. (|2j) to study the 
energy difference between two states that differ only by 
flux through loops of length i or greater 
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So, for example, if there is a state with zero flux, P n 
is maximized and it is the ground state. For a near- 
est neighbor model, this result proves a previously made 
conjecture[l3]. Interestingly, as a by-product of Eq. (|6]), 
we may also study the energetics of topological sectors 
that are characterized by loops that go all the way around 
the system (as discussed below). 

Enumeration of spin liquid states. In the quantum 
regime, the natural ground state may be a quantum spin 
liquid phase. Utilizing the above loop expansion, we can 
solve the energy minimization problem by first reducing 
the number of candidate ground states to a few spin liq- 
uid states and then comparing their energy. 

We begin by searching for spin liquid states among the 
set of translationally invariant Ising states with Q(ij\ = 
±1 within the unit cell, where (ij) are nearest neighbor 
bonds. That Qij is a pure phase results from the scale 
invariance introduced earlier. The restriction to Ising 
states then follows from the geometric frustration of the 
lattice and that Qij is a pair amplitude which breaks 
the U(l) gauge invariance down to ^[11]. We will not 
consider translationally non-invariant spin liquid states. 
This may be justified because the nature of the frustra- 
tion of the hyper-kagome lattice, being made of corner- 
sharing triangles, is similar to the kagome lattice where 
the Sp(7V) ground state is translationally invariant [ll) . 

The search for spin liquid states is then simplified by 
performing a symmetry analysis on the above set of mean 
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FIG. 3: The zero- and 7r-flux states of the hyper-kagome 
Sp(iV) quantum antiferromagnet. If an arrow points from 
site i to site j : then Qij = +1 and Qji = — 1. 



field states. By fixing the gauge, the symmetry of each 
state is then understood by studying its distribution of 
flux. The central result of this analysis [15| is that two 
Z2 spin liquid states and their 8-fold degenerate topolog- 
ical sectors are identified. The 8-fold degeneracy arises 
from the existence of tt flux through loops that go all the 
way around the system, and is expected when periodic 
boundary conditions are imposed on a three-dimensional 
Z2 spin liquid [16[. We find that these two states are dis- 
tinguished from each other by the existence or absence 
of 7r flux through the 10 site loop (the smallest even site 
loop on this lattice). We will therefore call them the 
zero-flux state and the 7r-flux state. 

The zero-flux state is the mean field ground state. This 
follows from the flux expulsion principle. In addition, we 
can show this explicitly by making use of Eq. [6] with 
P 5 = 18264 in the zero-flux state and P 5 = 18224 in the 
7r-flux state. The energy difference at small n is then 
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E Sp(N) [n-avLx}-E Sp{N) [0-fivLx} = — K 5 +O(K 6 ) > (7) 

The zero-flux state has zero flux through all of its smallest 
loops and is unique (up to topological degeneracy) in 
the translationally invariant sector. Note, the splitting 
between these states is small in the limit of small ft, a 
consequence of the long (length 10) loop that needs to 
be traversed to induce a splitting. 

Topological order. The topological order of the zero 
flux state, a central property of Z2 spin liquids, requires 
an eight fold degenerate ground state manifold with pe- 
riodic boundary conditions. In our construction, these 
eight degenerate states arise naturally. Consider, for 
example, the zero-flux state with an additional 7r flux 
through all z-axis "threads" (as highlighted in Fig. [T]) 
which wind around the system with L unit cells in the 
z-direction (L must be odd here). Since this state is oth- 
erwise the same as the state with zero flux through all 
threads, using Eq. (|6]) the energy splitting is 
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where (2L - 3)!!/(2L)!! « 1/ (2v^£ 3/2 ) for large L and, 
since there are only two loops of length 2L starting from 
the same site, P' r 



2L 



F 2L~ X = 4 /3. Note the exponen- 
tial dependence of AE on L is an explicit demonstration 
of the topological order of this phase [17j. Furthermore, 
since this is a three-dimensional Z2 spin liquid, a finite 
temperature transition to a paramagnetic state will occur 
in the inverted Ising universality class. This transition, 
however, is not in the spin sector and may not coincide 
with a cross-over in the spin structure factor from the 
high temperature CCP phase. 

Experimental signatures. The zero flux state has sev- 
eral features which distinguishes it from the CCP phase. 
Consider the energy integrated spin structure factor for 
the zero flux state shown in Fig. [2^ (this was calcu- 
lated following Ref. Illl ) and the static spin structure of 
the CCP phase (taken from Ref. l8|) shown in Fig|2j3. 
Notice that the spin liquid phase does not possess the 
long range dipolar spin correlations and its associated 
nodal structure along the [hhh] direction (and symmetry- 
related directions) — characteristic features of the CCP 
phase. The presence/absence of long range dipolar cor- 
relations clearly distinguishes these two phases. 

In the Z2 spin liquid phase, the elementary excitations 
are spin- 1/2 carrying spinions. Since the spin structure 
factor S(k, uj) measures spin-1 excitations, it vanishes un- 
less ou > e(q) + e(q — k) for any value of q where e(q) is 
the dispersion of a single spinon excitation. This defines 
the spinon- antispinon threshold 0. It is plotted as the 
bottom of the continuum given by the shaded (blue) re- 
gion in Fig. [2k. It should be entirely absent in the CCP 
phase. 

Semi- classical regime. Beyond the small k regime, the 
zero-flux state becomes unstable to magnetic ordering at 
k = n c = 0.4 (similarly the 7r-flux state becomes un- 
stable at k c — 0.8). This is remarkably large, given 
that magnetic order is expected to be more stable in 
three dimensions but that k c = 0.34 for a triangular lat- 
tice and n c = 0.53 for the kagome lattice [111]. The spin 
ordering pattern obtained upon spinon condensation in 
the zero-flux state is shown in Fig. [T^l. In contrast to 
state on the kagome lattice, this state has no 
local weather-vane modes[12]. Fluctuations within the 
classical spin manifold can only occur along an infinitely 
long thread with pattern BCBC. Thus, finite tempera- 
ture fluctuations may not select this state consistent with 
the spin nematic ordering and the absence of magnetic 
ordering found in numerics [8[. 

Conclusion. In this letter, we have presented tests for 
the existence of a quantum spin liquid phase in NIO. 
Using the large- TV Sp(iV) method, we proposed two can- 
didate ground states: a magnetically ordered state in the 
semi-classical regime and a Z2 quantum spin liquid state 
in a quantum regime. The magnetically ordered state 
has a q = coplaner spin ordering pattern which we ar- 
gue is unlikely to be entropically selected. Assuming NIO 



remains paramagnetic at low temperatures, it is still im- 
portant to find features distinguishing a quantum spin 
liquid from a classical cooperative paramagnet (CCP). 
Here we have discussed two such features for the Z2 spin 
liquid state. It should not have long ranged dipolar cor- 
relations and the associated nodal structures along the 
[hhh] direction (and symmetry related directions) that 
are the characteristics of the CCP phase. It should have 
a spinon-antispinon continuum with a threshold obeying 
a well defined dispersion. Future neutron scattering ex- 
periments on NIO at current temperatures and at lower 
temperatures, would be decisive tests for our predictions. 
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